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We investigate buckling of soft elastic capsules under negative pressure or for reduced capsule 
volume. Based on nonlinear shell theory and the assumption of a hyperelastic capsule membrane, 
shape equations for axisymmetric and initially spherical capsules are derived and solved numerically. 
A rich bifurcation behavior is found, which is presented in terms of bifurcation diagrams. The ener- 
getically preferred stable configuration is deduced from a least-energy principle both for prescribed 
volume and prescribed pressure. We find that buckled shapes are energetically favorable already 
at smaller negative pressures and larger critical volumes than predicted by the classical buckling 
instability. By preventing self-intersection for strongly reduced volume, we obtain a complete pic- 
ture of the buckling process and can follow the shape from the initial undeformed state through the 
buckling instability into the fully collapsed state. Interestingly, the sequences of bifurcations and 
stable capsule shapes differ for prescribed volume and prescribed pressure. In the buckled state, we 
find a relation between curvatures at the indentation rim and the bending modulus, which can be 
used to determine elastic moduli from experimental shape analysis. 

PACS numbers: 46.32. +x, 46.70.De, 62.20.de, 62.20.mq 



I. INTRODUCTION 

An elastic capsule consists of an elastic membrane en- 
closing a fluid phase. Elastic capsules are commonly met 
in nature, prominent examples are red blood cells or virus 
capsules. Artificial capsules can be fabricated by various 
methods [Ij, for example by interfacial polymerization at 
liquid droplets [2j or by multilayer deposition of poly- 
electrolytes [3], and have numerous applications. Sizes 
of capsules vary from the nanometer to the micrometer 
regime, and their mechanical properties depend on the 
fabrication process. For various applications, for exam- 
ple, if capsules are used as delivery and release systems, 
there is a need to characterize mechanical properties of 
capsules. The shape of capsules is approximately spher- 
ical but capsules are easily deformed by shear flow [4J, 
rotation [5], in adhesion [6l [7J, or by the application of 
local forces [F-TU]. Their deformation behavior also ex- 
hibits buckling instabilities upon decreasing the interior 
pressure or the enclosed volume [HI HHHH] or in adhesion 
[6]. All these deformation modes can potentially be used 
in experiments to infer material properties of capsules. 

Changes of the pressure inside the capsule by osmo- 
sis or mechanical means and changes in the capsule vol- 
ume represent the most basic deformation mechanisms 
for capsules with spherical rest shape. In this arti- 
cle, we study the collapse of a three-dimensional spher- 
ical capsule via the buckling instability into a fully col- 
lapsed state under negative pressure or for reduced cap- 
sule volume. The mechanical buckling instability sets in 
at the classical buckling pressure, which is well-known 
within linear shell theory for small displacements and an 
isotropic material Ll3j and has recently been extended to 
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shell materials with anisotropic shear response [T4]. The 
classical buckling pressure Pcb oc -E{Ho/Rof, where E 
is Young's modulus, Hq the initial membrane thickness 
and Rq its initial diameter, can be used in experiments 
to determine Young's modulus of the capsule material 

BUB. 

While the classical buckling pressure only marks the 
onset of the instability within linear shell theory and, 
thus, the limits of stability of a spherical capsule, it 
is much more difficult to calculate buckled shapes be- 
yond the critical buckling pressure because larger defor- 
mations require nonlinear theories and contact between 
originally opposite capsule sides has to be included in 
order to prevent self-intersections. Some nonlinear theo- 
ries have been applied to compute axisymmetric shapes 
with large deformations, for example in [15j under the 
assumption of isotropic tensions and for a more general 
case in [16j. Large buckling deformations have been con- 
sidered numerically based on triangulated surface models 
PI [TTHTQ] . In the framework of shell theory, however, a 
complete picture of the transition into buckled shapes as 
observed in the experiments is still lacking. 

In order to develop this picture for axisymmetric cap- 
sules, we use a nonlinear shell theory [20, 21j and assume 
hyperelastic capsule membranes. For hyperelastic ma- 
terials, a strain-energy function exists from which the 
tensions and bending moments can be deduced. It can 
be shown in general that solutions of the equations of 
force and moment equilibrium render the functional of 
total energy stationary [20j. This allows us to combine 
two tools, force equilibrium and principle of minimal en- 
ergy, in order to find different branches of stationary de- 
formed capsule shapes and then determine which capsule 
shape among different branches represents the global en- 
ergy minimum. At the classical buckling instability the 
branch corresponding to a spherical capsule loses stabil- 
ity and a bifurcation to buckled shapes takes place. This 
bifurcation is analyzed in detail beyond a linear stability 
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Figure 1. (Color online) Parameterization of the undeformed 
(left; always with index "0") and deformed midsurface 



analysis and in energy (and enthalpy) bifurcation dia- 
grams. We find that the buckling transition is discontin- 
uous in the energy diagrams and that the buckled shape 
becomes energetically favorable already at a smaller nega- 
tive pressure \pc\ < \pcb\ than the classical buckling pres- 
sure, where the spherical shape becomes mechanically 
unstable. 

We extend our approach to capsules with opposite 
sides in contact in order to prevent self-intersection at 
strongly reduced capsule volume. As far as we know, 
this has so far only been achieved for elastic rings [22j 
in two dimensions, but not for spherical shells in three 
dimensions. 

Furthermore, we analyze features of buckled shapes, in 
particular, the maximal curvature occurring at the edge 
of the indentation rim and find a relation between the 
maximal curvatures and the bending modulus, which can 
be used in experimental shape analysis. 



II. FINITE STRAIN SHELL THEORY 

A. Geometric Setup 

We consider axisymmetric capsules, which are under- 
going axisymmetric, torsionless deformations and whose 
membrane thickness is small compared to the other cap- 
sule dimensions. The midsurface of the undeformed 
capsule is parameterized by the curvilinear coordinates 
So G [0, Lq] and ^ [0, 27r) denoting the arc length 
measured along the meridian and the angle of revolu- 
tion, respectively. Its shape is determined by the func- 
tions ro(so) and zo{so) in cylindrical polar coordinates 
(see figure [l]) . 

In addition to the functions ro{so) and zo{so), the slope 
angle tpoi^o) is defined via the relations (see figure [T]) 



dro 
dso 



cos tpo and 



dso 



sin ipo . 



(1) 



Using this parameterization we calculate the principal 
curvatures of the midsurface. 



dipo , _ sin t/jQ 

tvgQ — ~, and tvfhQ — . 
dso ^ ro 



(2) 



Here, hZso denotes the curvature in the meridional direc- 
tion, and hZ(f)Q the curvature in the circumferential direc- 
tion. 

For the special case of capsules with spherical resting 
shape, the parameterization of the undeformed midsur- 
face is known analytically. 



^0(^0) = ^0 sin(5o/i^o) and 
^0(50) = i^o [1 - cos {so/Ro)] , 



(3) 



where the arc length sq ranges up to Lq = tt Rq. Accord- 
ingly, the slope angle and curvatures are given by 



i^o{so) = so/Ro and 



= V^o, (4) 



respectively. 

All quantities introduced so far in this section carry the 
index "0" because they refer to the undeformed capsule 
configuration. The midsurface of the deformed configura- 
tion is parameterized using analogous quantities without 
indices "0". Specifically, its shape is determined by the 
sought-after functions r{s), z{s) and the redundant tp{s). 
These functions have to satisfy the boundary conditions 
r(0) r{L) = z{0) = 0, V^(0) and ipiL) tt, corre- 
sponding to a closed capsule without kinks at its poles. 
The geometrical relations ([T]) and ^ can be transferred 
directly to the deformed midsurface by omitting all in- 
dices "0". 



B. Measures of Deformation 

Having defined the parameterization of the deformed 
and undeformed midsurface, we can now introduce mea- 
sures of deformation. Fibers oriented along the merid- 
ional and circumferential direction get stretched by the 
factors 



ds 
dso 



s\so) and A^^ 



r 

^0' 



(5) 



respectively. The function s{so) defined in this context 
describes the position s at which a particle can be found 
that was originally located at sq. In order to center the 
measures of stretch around zero, we define the meridional 
and circumferential strains = A5 — 1 and = A^ — 
1. The bending of the midsurface is measured by the 
meridional and circumferential bending strains 



K. = A. 



and = k.^ - k.^^ . (6) 



They are more suitable than a simple difference of de- 
formed and undeformed curvature because they lead to 
more simple constitutive equations, as we will see below. 



C. Elastic Law 

The strains measured by e^, e^. Kg and give rise 
to elastic tensions and bending moments in the capsule 
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membrane. Assuming that the capsule membrane con- 
sists of an hyperelastic material, there exists a surface 
energy density ws{es^ e^, i^s, K^j^)^ which measures the 
elastic energy that is stored in an infinitesimal patch of 
the membrane divided by the area that this patch takes 
in the undeformed configuration. In this paper, we will 
use a Hookean model [20j 



ws 



1 EHq 



{el^2ues 60 + 6^) 

^^Eb{K',^2uKsK^^KI) (7) 

with a (three-dimensional) Young modulus a bending 
modulus Eb^ and a Poisson ratio for a shell of (homo- 
geneous) thickness Hq. Note that the Poisson ratio of a 
two-dimensional membrane can take values — 1 < < 1, 
with V = 1 corresponding to an area-incompressible 
membrane. The product EHq is frequently called the 
surface Young modulus. In classical small strain theory 
of plates, the bending modulus can be expressed as 



Eb 



EH§ 



(8) 



12(1 

The bending energy contribution in the elastic energy ^ 
agrees to leading order in and (where Kg = Kq — i^sq 
and = i<i(j) — Hfj^^ ) and for an incompressible membrane 
with V = 1 with the commonly used Helfrich bending 

A2 



energy (/^g+z^^— cq) (see e.g. [181 [19]) with a spontaneous 
curvature cq = i^sq + • The Helfrich bending energy 
was originally proposed for two-dimensional liquids like 
vesicles, which differ qualitatively from the solid shells 
we are considering. 

It can be shown [20j that the meridional tension Ts 
and bending moment derive from the surface energy 
density via 



1 dws EHq 1 
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1 dws 
Xa, dKs 



1 - 
Eb 



(eg + ve^), 



(9) 
(10) 



Likewise, we obtain the circumferential tension r^j) and 
bending moment by analogous formulae with indices 
(j) and s interchanged. Tensions and bending moments 
are measured per unit length of the deformed capsule, 
which is the reason why prefactors 1/ A^ appear in these 
constitutive equations. 

Figure [2] shows on which faces of a membrane patch the 
tensions and bending moments act. Figure [2] also shows 
an additional transverse shear tension q which acts on 
the top (and bottom) side of the patch. It is constitu- 
tively undefined in our model because we did not incor- 
porate deformations in which the capsule's cross section 
gets sheared (i.e. fibers normal to the midsurface get ro- 
tated). However, the transverse shear tension is neces- 
sary to achieve force and moment equilibrium. Note that 
there does not act any transverse shear tension on the 
right side of the patch because of the assumption of ax- 
isymmetric, torsionless deformation. 




Figure 2. (Color online) Tensions and bending moments act- 
ing on the faces of a membrane patch 



D. Equilibrium Conditions 

Besides the geometric relations and constitutive equa- 
tions presented so far, conditions of force and moment 
equilibrium are needed to close the problem. These are 
three differential equations for tangential and normal 
force equilibrium and moment equilibrium, which take 
the form [20l M\ 



cosip 1 d(rrs) 

^ A 

r r as 

-p + r0 + his Ts + 

co^ijj 1 d(r rris 
- 



1 d{rq) 
r ds ' 



ds 



(11) 
(12) 

(13) 



in the absence of external tangential force and torque 
densities and for a constant pressure p inside the capsule. 



III. PRINCIPLE OF STATIONARY ENERGY 

A. Variation of Energy Functionals 

Another approach to find stable configurations of cap- 
sules for fixed but altered volume V ^Vq is to minimize 
the functional of stored elastic energy F = J ws dA^ by 
calculus of variations. The constraint of fixed capsule vol- 
ume V is handled by introducing a Lagrange multiplier p 
and extremizing the enthalpy G = F — pV instead. The 
principle of stationary total energy [20 states that the 



solutions of the equilibrium conditions (11) to (13) with 



given pressure p render the enthalpy G stationary. 

In the case at hand, this can be verified by using 
the standard procedure of calculus of variations. Using 
dAo = 2 TT ro dso for the area element of the undeformed 
midsurface and V = Jnr'^ dz = /ttt^ sim/j ds = 
/ irr^ Xs sin ip dso as an integral expression for the cap- 
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Figure 3. (Color online) Ansatz for the configuration with 
opposite sides in contact. Along the contour, the important 
coordinates sq = Lc and sq — Lq/2 are marked. 



sule volume, we arrive at the enthalpy 



G 



(27TroWs—p7Tr'^XsSimjj^ dso, (14) 



for which we have to find stationary points. The Euler- 
Lagrange equations are derived in appendix [A] and turn 
out to be 



cos^^ 1 dirrris) 
^0 ~ r — JZ ^ 



= 



r 

1 d{rrs 



r ds 
cos iIj 



ds 



1 



r 



Ks q 



-Ts tanV^+ - p- 

2 cos 



(15) 
(16) 
(17) 



They coincide with the general equations of tangential 
force and moment equilibrium (11) and (13), but the 



normal force equilibrium equation (12) is replaced by the 



algebraic expression (17) for q. Instead of three differen- 



tial equations in the general case, we now end up with 
only two differential equations and one algebraic relation. 
This means we have found a first integral of the general 
equations of equilibrium (which is only valid for capsules 
with spherical resting shape under uniform pressure). In- 
serting the algebraic expression (17) into the differential 



equation (12) for g confirms that the solution found here 



satisfies the general equations of force equilibrium. 



B. Shapes with Opposite Sides in Contact 

Among the solutions for negative pressure, unphysi- 
cal solutions that exhibit self-intersection can occur. For 
two-dimensional shells with circular resting shape (cylin- 
der, ring), Flaherty et al. addressed this problem in [22j. 
They found that the opposite sides of an elastic ring 
touch at one point when the negative pressure reaches 
a certain threshold Pci- Lowering the pressure further, 
the curvature at the point of contact decreases until it fi- 
nally becomes zero at a second critical pressure Pc2- For 
pressures lower than Pc2, the contact area is a straight- 
line segment, which increases in length with decreasing 
pressure. 

In order to generalize these results to three- 
dimensional spherical shells, we make the ansatz of a 
top/down symmetric deformed configuration with flat 



circular areas around the poles in contact with each 
other, see figure [sj Considering only top/down symmet- 
ric solutions, it is sufficient to set up the energy func- 
tional for the lower hemisphere only. The flat part with 

< 5o < I/c? is referred to as part I, and the remaining 
part with Lc < sq < as part II. In order to treat 
variations with respect to the boundary Lc between parts 

1 and II properly, the enthalpy functional is split up into 
two corresponding parts 

G = W -pV - ii/^z 

dso{27rroWs) (18) 

rLo/2 

+ / dso (27rro i^^s- — pTrr^ sin?/; As — /i As sin?/;), 

where the volume constraint and the condition = Az = 
J dz = J As sin?/; d^o are already incorporated via the 
Lagrange multipliers p and /i. 

Performing the first variation involves the following 
problems: In part I, the function r{so) must be var- 
ied, and in part II the functions r(so) and ?/^(5o). Addi- 
tionally, the integral boundary Lc must be varied, which 
yields continuity conditions at sq = Lc. The results can 
be summarized as follows. 

In part I, we obtain one Euler-Lagrange equation 



djrTs) 
ds 



for So e [0, Lc]. 



(19) 



In part II we obtain a system of two Euler-Lagrange equa- 
tions 



: 



Q 



cos ip 

ms - 

r 

1 djrrs) 
r ds 

—Ts tan if) - 



1 dirrris) 
r ds 



r 



1 



r 



oP T 

2 cos ip 



1 



27r r cosip 



(20) 
(21) 
(22) 



for So G [Lc, Lo/2]. The Euler-Lagrange equations for 
part II are very similar to those obtained above for non- 



intersecting shapes, except for the term 



in the 



27r r cos ip 

algebraic expression for q{so). Nevertheless, it can be 
easily shown that relation ( [22| ) also satisfies the general 
equation (12), although it differs from the previous solu- 



tion by the additional term with the Lagrange multiplier 
/i. 

In order to fit the solutions of part I and II together, 
continuity conditions must be derived, which arise from 
the variation with respect to the integral boundary Lc- It 
turns out that most quantities of interest are continuous 
at So = Lc, namely r, ?/^, A^, z^^, r^, rrii where an 
index i stands for s or cj). Only for the transverse shear 
tension q no continuity condition derives because of the 
lack of a constitutive relation between the shear tension q 
and elastic deformations in the present formulation. We 
just know that the value q{Lc) depends via the algebraic 
relation (22) on the Lagrange multiplier /i. 
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IV. SHAPE EQUATIONS 

The shape of the deformed capsule is governed by the 
equihbrium conditions, constitutive equations and geo- 
metric relations presented above. They can be rear- 
ranged to form a system of first order differential equa- 
tions, known as the shape equations^ 



r'iso) 


= As costjj 


z'iso) 


— As sin i/j 


i^'iso) 


— As l^s 


<iso) 


= As I cos — — — + 
V r 


m;(so) 


— As cos w — 


V r 


q'iso) 





(23) 



The first three equations in this system follow from the 
geometric relations ([T]) and Q (without index "0") and 
the last three from the equilibrium conditions (11) to 



(13). The change of variables from s to sq was accom- 
plished by the relation ds = As d^o, see ([5|. Usage of the 
algebraic expression (17) would reduce this system by 



one equation, but turns out to be numerically impracti- 
cal because of singularities when approaches 7r/2. 

In order to close the above system of shape equations, 
all functions appearing on the right hand side must be 
expressed in terms of the basic functions r, ip, Tq^ rris 
and q. Exploiting geometrical relations, the definitions 
of the strains and the constitutive equations, we find the 
set of relations 

H , As = (1 - ) A0 — -- - iy{X^ - 1) + 1, 

^ H/ rLc 



r 







>^cb = — , r0 = _ ^((A0 - 1) + i^(As - 1)), 



sin ip — sin tpo 



Ks = -^Xcpms-i^K^ 
Eb 



m^ = Eb {K^ + u Ks) , K.s ^ 
As As 



(24) 



which close the system (23) 



After changing variables from s to 5o, the boundary 
conditions for a closed capsule without kinks at its poles 
are r(0) = r(Lo) z(0) ?/;(0) and MLq) tt, 
while the spherical reference shape satisfies ([3| and (|4]). 
These boundary conditions cannot be directly used to 
determine boundary values for all strains and curva- 
tures using ( [24| ) because some of the resulting expres- 
sions are ill-defined as ratio of two vanishing quantities. 
These boundary values can be evaluated analytically us- 
ing L 'Hospital's rule and symmetry arguments (physi- 
cally significant functions should be either odd or even 
along the extended contour sq G [—Lq^ Lq\) in the limits 
5o ^ and sq Lq. Using this procedure boundary 



values at 5o = and sq = Lq can be written as 



Ae — Xrh — 



^i^o-rs(l 



1 



EBil^iy) RoXs 



0. 



(25) 

Now, the system ([23]) can be solved numerically. In 
order to introduce dimensionless quantities, we choose Rq 
as the length unit and E Hq as the tension unit. A simple 
shooting method and a multiple shooting method [23l 
l24j . both with parameter tracing, were used to compute 
capsule shapes for progressively lowered pressure p < 0. 

For top/down symmetric configurations it is sufficient 
to integrate from the south pole 5o = to the equa- 
tor So = Lq/2. In this case, the boundary conditions 
r(0) = 0, z{0) = 0, V^(O) = 0, q{0) = 0, ^{Lo/2) = 7r/2 
and q{Lo/2) = must be satisfied. Whereas the condi- 
tions concerning r, z and ip are obvious from geometry, 
the condition q{0) =0 holds because of the algebraic re- 
lation (17), and q{Lo/2) = because of the symmetry of 
the configuration: If q{Lo/2) ^ 0, the lower hemisphere 
would pull the upper one to the inside (or outside, de- 
pending on the sign of g(I/o/2)), which contradicts the 
top/down symmetry. 

For solutions that are not top/down symmetric, we use 
shooting to a fitting point [23j because it is numerically 
impractical to integrate into a removable singularity. The 
appropriate boundary conditions are r(0) = 0, z(^) = 0, 
V^(O) = 0, ^(0) = 0, r(Lo) = 0, V^(Lo) = ix and g(Lo) = 
in this case. The apparent problem that there are seven 
boundary conditions to a system of six first order differ- 
ential equations is resolved by the algebraic relation (17). 



It effectively renders one of the boundary conditions to 
q obsolete (in other words, the differential equation for q 
with both boundary conditions concerning q is obsolete, 
which leads to a system of five differential equations with 
five boundary conditions). 

In the case of top/down symmetric configurations with 
opposite sides in contact, the solutions of part I and II 
must be fitted together according to the continuity condi- 
tions. In part I, the only degree of freedom is the function 
r(5o) since 2;(5o) = '^(sq) = because of the geometrical 
restrictions. The shape equations for part I can be ob- 
tained from the equilibrium condition (19), together with 
constitutive and geometric relations. 



/(5o) =As 



A. 



Ro 



(26) 



For reasons of symmetry, As(so) must be an even function 
along the extended contour sq G {—Lq, Lq). Therefore, 
Xs{0) = 0. The boundary condition for this part of the 
shape equations is r(0) = 0. The stretch at the pole, 
As(0), is free and serves as a shooting parameter. 
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After a numerical solution of p6| ) has been found, 
the Hookean constitutive equations can be used to cal- 
culate the tension Tc = Ts{Lc) and bending moment 
rric = ms{Lc), which serve together with Vc = r{Lc) as 
initial values for the integration on part II. The shape 
equations for part II can be adopted from ([23|. Only 
the boundary conditions must be changed to r(Xc) = ^c? 
z{Lc) 0, tp{Lc) 0, TsiLc) Tc, ms{Lc) = rric, 
z{Lo/2) = 0, V^(Lo/2) = 7r/2 and q{Lo/2) = 0. The 
shooting parameters to satisfy the three conditions at 
the far end are A5(0), q{Lc) and Lc. 

In the case of top/down asymmetric configurations 
with opposite sides in contact, we assume the region in 
contact to be a spherical cap with radius R. A region 
< So < I/c around the south pole has to change its 
curvature from 1/Ro to —1/R (note the sign change) to 
match a region Lq — Lc2 ^ sq < Lq around the north pole. 
We allow the region around the south pole to contract in 
the meridional direction by the factor Xs{so) = A^^ (con- 
stant on < 5o < i^c)- The region around the north pole 
is allowed to contract by the factor Xs{so) = (con- 
stant on Lq — Lc2 < ^0 < Lq) which may be different 
from Xd. Since the deformed configurations of these two 
regions have to match exactly, we have the constraint 
Lc2 = Lc Xd/Xt. 

Thus, the shape of the contact region of the deformed 
capsule is determined by the four parameters R^ Lc^ A5 
and At. The shape of the non-contacting part of the 
capsule is again described by the shape equations (23). 



Since it is a system of six equations, and we have four 

additional parameters that are not known a priori, we are 

able to satisfy ten continuity conditions in total. That is 

just enough for the most important quantities r, 

and rris at both ends Lc and Lq — Lc2 of the integration 

interval. 

The explicit boundary conditions for the shape equa- 
tions ([23| are given by 

r{Lc) = R sm{Lc Xd/R) 
z{Lc) = R{cos{LcXd/R)-l) 
^{Lc) = -LcXd/R 
r(L,) = ^^(A,-l + KA0-l)) 
ms{Lc) = EB^{-Xd/R - 1/Ro + ^{-X^/R - l/R^)) 
with X^ = R sm{LcXd/R)/{Ro sm{Lc/Ro)) 

(27) 

for the starting point and 



r{Lo- 


Lc2) 


z{Lo - 


Lc2) 


i>{Lo- 


Lc2) 


t{Lo- 


Lc2) 


ms{Lo - 


Lc2) 


with 


X(h 



EHo 1 



for the end point of integration. 



(28) 



Note that taking A5 to be constant in the regions in 
contact is a strong simplification. A more involved the- 



ory would incorporate equations similar to ( 26 ) to deter- 
mine the in-plane-displacements. However, we observed 
that the solution branches produced with this simplified 
method fit neatly in the bifurcation diagrams. 



V. BIFURCATION BEHAVIOR 

As the surface Young modulus E Hq serves as the ten- 
sion unit, there are only two elastic parameters left to 
vary, the Poisson ratio v and the dimensionless bending 
modulus 



Eb 



Eb 



RlEHo 12{l-u^)Rl' 



(29) 



where we used (|8| for the bending modulus of elastic 
plates. We will present bifurcation diagrams for a Poisson 
ratio of u = 0.5 and dimensionless bending moduli of 
Eb =0.01 and 0.001. These values correspond to relative 
shell thicknesses of Hq/Rq = 0.3 and 0.095, respectively. 



A. Spherical Solution Branch 

The trivial branch of spherical solutions (branch 1 in 
the bifurcation diagrams below) can be calculated ana- 
lytically because of its high symmetry. Solving the equi- 
librium equations (Ht]) to dlsl) for a sphere with radius R, 



it is straightforward to show that these equations reduce 
to the modified Laplace- Young equation 



p = 2 tzr^ 



(30) 



with K = 1/R the isotropic curvature and r = Ts = T(j) 
the isotropic tension. The Laplace- Young equation deter- 
mines the new radius of the capsule for given pressure. 
Using the Hookean constitutive relation to express r in 
terms of the isotropic stretch A5 = A^ = i^/i^o? this con- 
dition can be rewritten as 



2 EHo 



{R - Ro) 



0. 



(31) 



The solution to this quadratic equation is the radius- 
pressure relation of the spherical branch (called branch 1 
below) 



Ri 



EHn 



p(l-zy) 



± 



EHn 



2 EH, 







Ro. (32) 



where the + branch holds for p < and the — branch 
for p > 0, which can easily be inferred from requiring 
limp^o R = Rq- This radius-pressure relation determines 
the deformed configuration completely, and all physical 
properties, like volume, strains, tensions and stored elas- 
tic energy, can be calculated in turn. 
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In the following we focus on buckling shapes for neg- 
ative pressure p < 0. For positive pressure p > 0, where 
capsules are stretched, the spherical branch 1 represents 
the only equilibrium shape. 

B. Stability Criteria 

We have shown that the equations of force and moment 
equilibrium render the functional of total energy station- 
ary When the shape equations are solved for negative 
pressure, various solution branches with reduced capsule 
volume V < Vq can be found, which can represent local 
minima or maxima of the energy functional. We display 
the total energy or enthalpy of each solution branch in an 
energy bifurcation diagram as a function of the capsule 
volume V or the pressure p, respectively. The principle of 
minimal total energy allows us to determine the globally 
stable branch as the branch of minimal energy among 
all stationary shapes. Shape transitions such as buckling 
occur where two branches intersect. 

If the capsule volume is given, the stored elastic energy 

F = 271 tqWs dsQ (33) 

must be minimal. This criterion has experimental signif- 
icance if the capsule volume cannot change because the 
encapsulated liquid is incompressible and the membrane 
is impermeable. In cases of a semipermeable capsule 
membrane, it is also reasonable to consider the volume 
fixed because the relaxation into the equilibrium shape 
happens on much shorter time scales than the diffusion 
of the inner liquid through the membrane. In the corre- 
sponding bifurcation diagram we display the stored elas- 
tic energy F as a function of the volume V . 

On the other hand, if the capsule is filled and sur- 
rounded by gases, the pressure difference p is prescribed 
rather than the capsule volume. In this case, configura- 
tions with minimal enthalpy 

G = F — pV = / {27r To w s — p TT r"^ Xs sinilj) dso 
Jo 

(34) 

are energetically preferable. In the corresponding bifur- 
cation diagram we display the elastic enthalpy G as a 
function of the pressure p. Note that G{p) is the Legendre 
transform of F{V), since the relation p = dF/dV holds 
(which was verified numerically for all solution branches 
presented here). Solution branches that lie lowest in the 
F{V) bifurcation diagrams need not necessarily coincide 
with the lowest branches in the G{p) diagrams. 

Finally, it is also useful to analyze the relation between 
pressure p and volume V for stationary capsule shapes. 
Branches that exhibit an unusual pressure-volume rela- 
tion with dp/dV < are inherently unstable if the pres- 
sure is given instead of the volume [25j. To see that, we 
consider a water filled capsule connected to a reservoir 
of water. The pressure in this system can be prescribed. 



Now, if we try to change the capsule volume by lowering 
the pressure by an amount dp < 0, the capsule grows by 
an amount dV > 0, i.e. water flows from the reservoir into 
the capsule. The loss of water in the reservoir typically 
leads to an even lower pressure and thus, the equilibrium 
is unstable with repsect to volume changes. This insta- 
bility is also reflected by a negative second derivative of 
the free energy, dPF/dV'^ = dp/dV < and a horizontal 
p{V) curve with dp/dV = marks the onset of such an 
instability. 

Therefore, we have three criteria of stability for capsule 
shapes: 

(i) minimal energy F for fixed capsule volume V ^ 

(ii) minimal enthalpy G for fixed pressure p, and 

(iii) a monotonously decreasing pressure-volume rela- 
tion dp/dV < is sufficient for an unstable shape 
for fixed pressure. Thus, the criterion dp/dV > 
is only a necessary condition for stability: Config- 
urations with dp/dV > can still be unstable with 
respect to deformation modes which do not change 
the volume. 

In fact, criterion (iii) can be generalized by using a gen- 
eral bifurcation theorem that has been proven in Ref. [26j . 
This allows us to make a further statement about the in- 
stability of shapes beyond points where a monotonously 
decreasing p(y) curve becomes vertical: If the p{V) curve 
at such a point is open to the left, i.e., the following 
lower part of the p{V) curve has again a positive slope 
dp/dV > 0, this lower part must also be unstable. Fur- 
thermore, this is an instability with respect to a volume- 
preserving mode. This generalization also demonstrates 
that a positive slope dp/dV > is not sufficient for sta- 
bility. 



C. Bifurcation Diagrams 

With these criteria of stability, we analyze the bifurca- 
tion behavior of elastic capsules by investigating the dif- 
ferent branches in three types of bifurcation diagrams: 1) 
In the F{V) diagram we study buckling by reducing the 
capsule volume by using criterion (i). 2) In the p{V) di- 
agram we can identify unstable shapes as monotonously 
decreasing branches dp/dV < according to criterion 
(iii). 3) In the G{p) diagram the Legendre transform 
G{p) = F{V{p)) — pV{p) allows us to investigate buck- 
ling under negative pressure according to criterion (ii). 

The nomenclature of solution branches is summarized 
in table |Tj In particular, we will compare our results to 
classical buckling theory fi3\ \20\ |27j, which predicts a 
critical negative buckling pressure 

= -^ n2 /fo^° 2 T = -^VeH^/RI (35) 
Ri a/12(1 - u^) 
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Figure 4. (Color online) Bifurcation diagram for given volume of a capsule with Eb — 0.001. On the dashed lines of branch 2 
and 3, the capsule intersects itself. Branch 3' continues winding up in the diagram, as indicated by the "etc 



P/Pcb 




Figure 5. (Color online) Pressure-volume relation of a capsule with Eb = 0.001 
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Figure 6. (Color online) Bifurcation diagram for given pressure of a capsule with Eb = 0.001. 



Branch 


Configuration 


Example 


1 


spherical 


o 


2 


simply buckled (asymmetric) 


Q 


2' 


asymmetrically crumpled 


O 


2c 


simply buckled with contact 


/~\ 


3 


symmetrically buckled 


CD 


3', 3" 


symmetrically crumpled 


OO 


3c, 4c 


symmetrically buckled with contact 


f> 0,0 o 



Table I. Denomination of different branches in bifurcation 
diagrams. Branches are classified according to the num- 
ber of bulges as buckled (one or two bulges) or crumpled 
(more than two bulges, names with primes) and according to 
top/down mirror-symmetry (shapes 1,3,4 are top-down sym- 
metric, shapes 2 are asymmetric). Shapes with subscript c 
exhibit contact of originally opposite sides; shape 3c develops 
an additional dimple at the sides as opposed to shape 4c (for 
Eb = 0.001). 



for a spherical capsule. A corresponding critical volume 
can be obtained as 



Vr 



47r 



cb 



(36) 



by using pcb in the radius-pressure relation (32) of the 



spherical branch 1. 

We start with the bifurcation behavior of thin capsules 
with Eb = 0.001, as shown in the diagrams [i] to [i] The 
F(V) diagram (figure |4| reveals that the simply buck- 
led configurations of branch 2 are energetically favorable 
for volumes V < Vc, where Vc denotes a critical volume 
which is Vc ~ 0.944 Vb in this case. At the critical vol- 
ume Vc^ the spherical branch 1 and branch 2 of simply 
buckled configurations intersect. Note that this volume is 
larger than the critical volume Vcb within classical buck- 
ling theory, which is Vcb ^ 0.914 Vq for this case. A shape 
transition between branches 1 and 2 at Vc is discontin- 
uous and involves an energy barrier. The upper part of 
branch 2 most likely represents the unstable transition 
state ai V = Vc between a spherical shape and the stable 
lower part of branch 2. Therefore, the energy barrier can 
be estimated by the energy difference between the upper 
and lower parts of branch 2 dit V = Vc- More detailed 
stability considerations are given in section [V D| below. 

For volumes V < 0.14^0 7 these configurations start 
to intersect themselves (dashed line). Simultaneously, 
a branch 2c starts to exist with simply buckled config- 
urations with opposite sides in contact. Although the 
method used to obtain this branch incorporated some 
simplifications, branches 2 and 2c connect neatly in the 
diagrams. In this domain, the ansatz with opposite 
sides in contact produces solutions with higher energies 
(branch 2^ compared to the dashed part of 2), i.e. the 
self-intersecting solution branch raises in the F{V) dia- 
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gram when it is forced to satisfy the physical constraints. 

For volumes V < Vc^ the energetically next best con- 
figurations according to (i) are the symmetrically buck- 
led ones of branch 3. They exhibit self-intersection for 
volumes V < O.^IVq. At higher energies, crumpled con- 
figurations can be found in branch 3' that is connected 
to branch 3. They correspond to local minima, saddle 
points or maxima of the energy functional. Several more 
crumpled shapes can be found in this domain of the bi- 
furcation diagram, but they are not shown for the sake 
of clarity. Similar crumpled configurations have been ob- 
served for small volumes in simulations using triangu- 
lated surface models, in particular at high compression 
rates [17J. Trapping in metastable crumpled shapes can 
contribute to this behavior. 

For given pressure, the p{V) diagram (figure [5| and 
G(p) diagram (figure [6| can be analyzed in order to 
identify unstable solutions according to criteria (ii) and 
(iii). The simply and symmetrically buckled solutions of 
branches 2 and 3 (without contact of opposite sides) ex- 
hibit a negative slope in the p{V) diagram. This means 
that they are mechanically unstable with respect to vol- 
ume reducing deformations for given pressure according 
to criterium (iii) , although they are most stable for fixed 
volume in the F{V) diagram. According to the general- 
ization of criterion (iii) based on Ref. [26j also the lower 
part of branches 2 and 3, which lie beyond the turning 
point and have a positive slope in the pjV) diagram at 
volumes V > Vcb, (see also inset in figure pi) are unstable, 
however, with respect to volume-preserving modes. The 
energy diagram G{p) confirms this result; the unstable 
branches 2 and 3 lie above the trivial solution branch 
1. On the other hand, the buckled configurations 2c and 
3c with opposite sides in contact are mechanically stable 
again, with dp/dV > 0, and 2c is energetically preferable. 

At the critical pressure Pc ^ —0.23 \pcb\^ the spherical 
branch 1 and branch 2c of simply buckled configurations 
with opposite sides in contact intersect in the G{p) dia- 
gram, and shapes 2c become energetically favorable for 
given pressure. 

The negative critical pressure \pc\ is much smaller 
than the classical buckling pressure \pc b \ , which is Pcb ^ 
— 0.126 £^i^o/^o according to eq. ( [35| ). However, the 
classical buckling pressure fits very well to the point, 
where branches 2 and 3 emerge from the spherical so- 
lution branch in the p{V) diagram (see inset in figure [sj 
where it is indicated by a horizontal line). The classical 
buckling pressure Pcb is the pressure where the unbuckled 
spherical configuration 1 becomes unstable with respect 
to buckling [13, 14J, and a spontaneous transition to the 
unstable simply buckled branch 2 occurs. 

At the critical pressure Pc, on the other hand, shape 
2c becomes energetically favorable as compared to the 
spherical branch 1. A shape transition between both 
branches at Pc is discontinuous and involves an energy 
barrier. The upper unstable branch 2 of a simply buck- 
led shape without contact most likely represents the tran- 
sition state between a spherical shape and branch 2c at 



p = Pc Therefore, the energy barrier can be estimated by 
the energy difference between the upper unstable branch 
2 and the lower stable branch 2c at p = Pc- 

The bifurcation diagrams for a thick capsule with 
Eb = 0.01 (figures [t] to |9| look qualitatively similar. 
Again, the simply and symmetrically buckled solution 
branches 2 and 3, respectively, are energetically prefer- 
able for given volume, but exhibit a negative slope in the 
p{V) diagram (figure [8|, and are unstable and energeti- 
cally unfavorable for given pressure. 

For given volume the F(V) diagram (figure[7|) shows 
that the simply buckled configurations of branch 2 are 
energetically lower than spherical shapes for V < Vc with 
a critical volume Vc ^ 0.79Vo, which is again larger than 
the critical volume Vcb ^ 0.77^0 from classical buckling 
theory. 

For given pressure p, the G{p) diagram (figure[9| shows 
that simply buckled configurations 2c with opposite sides 
in contact are energetically favorable as compared to 
spherical shapes of branch 1 for p < p^ where the critical 
pressure is Pc ^ —0.46 \pcb\' Also for Eb = 0.01, the neg- 
ative critical pressure \pc\ is much smaller than the clas- 
sical buckling pressure \pcb\^ which is Pcb ^ —OAEHq/Rq 
in this case and the classical buckling pressure fits to the 
point, where branch 2 emerges from the trivial spherical 
solution branch (see inset in figure [sj where it is indicated 
by a horizontal line). 

For = 0.01, there is a visible gap between the buck- 
led branches 2 and 3 and their respective continuations 
2c and 3c with opposite sides in contact. This gap was 
already present in the diagrams for Eb = 0.001, but 
much smaller. It is assumed to be closed by configura- 
tions with point contact of north and south pole. In the 
case of branch 3, an analogous behavior like that of elas- 
tic rings is expected [22j. The curvature at the point of 
contact is expected to decrease, until it finally becomes 
zero and hence fulfills the continuity conditions for circu- 
lar areas in contact. However, the shape equations when 
point contact of north and south pole is enforced turn out 
to be hard to solve numerically, because the transverse 
shear tension diverges at the poles. 

At higher energies, there are again configurations with 
several bulges (branches 3' and 3"). In the F{V) diagram 
(figure [t]), they lie lower than the trivial solution branch. 
In contrast to the results for the capsule with Eb = 0.001, 
branch 3" does not have multiple turning points and is 
not connected to the continuation 3' of the symmetrically 
buckled branch within the scope of our diagrams. How- 
ever, branches 3' and 3" might join at higher energies 
and lower pressures. Notably, these solution branches lie 
lower than branch 1 in the G{p) diagram (figure|9|, which 
is a qualitative difference to the results for = 0.001. 

The F{V) and p{V) diagrams are in good agreement 
with previous work based on triangulated surface models: 
In Ref. [19j, ap(V) relation has been obtained, which also 
shows a uniform shrinkage of the capsule (our branch 1) 
for small volume reduction followed by a jump into an ax- 
isymmetric simply buckled configuration (our branch 2) 
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Figure 7. (Color online) Bifurcation diagram for given volume of a capsule with Eb = 0.01. The dashed lines in the gaps 
between 2/2c and 3/3c could not be calculated numerically, but we expect some configurations to exist in these domains. 




Figure 8. (Color online) Pressure- volume relation of a capsule with Eb = 0.01. The dashed lines in the gaps between 2/2c and 
3/3c could not be calculated numerically, but we expect some configurations to exist in these domains. 
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Figure 9. (Color online) Bifurcation diagram for given pressure of a capsule with Eb = 0.01 



with the same p{V) behavior as branch 2. Furthermore, 
Ref. fT8] contains an F(y) diagram in which the spher- 
ical and simply buckled solution branches are shown. 
They reveal qualitatively the same F{V) relation as our 
branches 1 and 2. 



D. Buckling Bifurcation 

The key features of the bifurcation diagrams presented 
above are drawn schematically in figure [To) They allow to 
construct a complete picture of the bifurcation behavior 
at the buckling transition. 



On the left of figure 10, two different bifurcation sce- 
narios in the domain where branch 2 emerges from 1 are 
drawn schematically. We see that the part of branch 
2 with an inward dimpled south pole emerges continu- 
ously from the spherical branch 1 and first runs to the 
right, i.e. to higher capsule volumes. After a turning 
point, the branch runs to the left, i.e. to lower capsule 
volumes, while the south pole buckles more and more 
inwards. The upper part of branch 2 up to the turn- 
ing point (shaded dark red) lies at higher energies than 
the spherical branch and consists of unstable stationary 
shapes. This can be seen from the p{V) diagram where it 
corresponds to the lower part of branch 2 beyond the ver- 
tical turning point and which are unstable with respect 
to a volume-preserving deformation mode according to 
the generalization of criterion (iii) [26j. 



Branch 2 intersects the spherical branch 1 at the criti- 
cal volume V = Vc. For V <Vc the lower part of branch 
2 is the energetically preferable solution branch. A shape 
transition between branches 1 and 2 at Vc is discontinu- 
ous. If the capsule wants to switch from the metastable 
branch 1 to 2 for V < Vc (vertical arrow) an energy bar- 
rier must be overcome. The upper part of branch 2 rep- 
resents the unstable transition states between a spherical 
shape and the stable lower part of branch 2. Therefore, 
the energy barrier can be estimated by the energy differ- 
ence between the upper and lower parts of branch 2. 

The behavior when branch 2' emerges from the spher- 
ical branch 1 at the volume V = Vcb is quite different. 
It emerges with egg-like configurations continuously from 
the trivial branch and runs directly to the left. If the cap- 
sule passes this point in the bifurcation diagram during a 
progressive reduction of its volume along the metastable 
branch 1, it is allowed to switch from branch 1 to 2' con- 
tinuously. Thus, there is no energy barrier to be overcome 
in this scenario, and the trivial branch is supposed to be 
unstable. 

Details of a realistic buckling process for given pres- 
sure can be constructed from the p{V) and G{p) dia- 
grams. The middle and right parts of figure [To] schemat- 
ically show the key features concerning the simply buck- 
led states 2 and 2c. The decreasing part of branch 2 
with dp/dV < is mechanically unstable with respect to 
volume reduction (shaded light red). According to the 
generalization of criterion (iii) based on Ref. [26j also the 
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Figure 10. (Color online) Left: Schematic drawing of two different bifurcation phenomena (indicated by the arrows) for given 
volume V < Vc. Middle and right: Schematic drawing of a realistic buckling process for given pressure. The critical pressure 
Pc is defined by the crossing of branches 1 and 2c in the G{p) diagram. As a consequence of the Legendre transformation, the 
two regions shaded gray in the p{V) diagram have the same area at p = pc (Maxwell construction). 



lower part of branch 2 (shaded dark red), which has a pos- 
itive slope for a small volume range V > Vcb in the p{V) 
diagram, is unstable with respect to a volume-preserving 
deformation mode. It corresponds to the unstable upper 
branch 2 in the F(y) diagram. In the G{p) diagram both 
corresponding unstable parts of branch 2 join to give the 
energetically unfavorable upper branch. 

Therefore, the spherical shapes of branch 1 become 
mechanically unstable at the classical buckling pressure 
Pcb, where the unstable branch 2 and branch 1 merge 
in the p{V) diagram. Then, a small dimple caused by 
fluctuations or agitations can grow spontaneously and 
the capsule finally ends up in a fully collapsed stable 
configuration 2c at the given pressure running along the 
dashed path in the p{V) diagram. The same process is 
indicated in a schematic G{p) diagram on the right of 
figure [To] as a dashed line. 

Shape 2c becomes energetically preferable already at 
a much smaller negative pressure Pc, i.e., \pc\ < \Pcb\, 
where branches 1 and 2c intersect in the G{p) diagram. 
A shape transition between branches 1 and 2c at Pc is dis- 
continuous: Changing onto branch 2c at p = Pc does not 
correspond to a spontaneous snap-through into a fully 
buckled shape 2c because shape 1 remains mechanically 
(meta)stable in that region, because dp/dV > and no 
other branches are merging/intersecting in the p{V) dia- 
gram at Pc- A finite dimple has to form by fluctuations 
or agitations to induce buckling, and this is associated 
with an energy barrier. The unstable transition state of 
this process can be the unstable shape 2 at the same pres- 
sure. This process is shown as a solid path in the p{V) 
diagram. We note that, as a consequence of the condition 
of equal enthalpies G at p = Pc, this solid path can be 
obtained by a Maxwell construction: At p = pc, the two 
regions shaded gray in the p{V) diagram have the same 
area. 

We conclude this section by providing estimates for 
buckling pressures for some synthetic and biological cap- 
sules. For synthetic capsules made from typical soft 
materials we expect a Young's modulus in the range 
- 100 - 1000 MPa, thicknesses i^o 10 - 50 nm, 



and micrometer sizes Rq ~ 500 nm, see for example Refs. 
|8l [To] for different synthetic capsules. This results in 
typical classical buckling pressures \pcb\ ~ 0.1 — 1 GPa 
(in accordance with measurements in [8j). Such materi- 
als have a small dimensionless bending modulus Eb ^ 
0.00005 — 0.001 corresponding to thin shells, see figures 
|4)[5l and[6l 

Many biological materials, such as virus capsids have 
very similar material characteristics but can be smaller: 
In Ref. [14] E 1 GPa, Hq 2nm, and 10 nm 

has been used for virus capsids, which gives similar di- 
mensionless bending modulus Eb ~ 0.004 and a similar 
buckling threshold \pcb\ = 0.5 GPa. 

Somewhat different are soft biological capsules such 
as red blood cells with a shell made from lipid bilayers, 
which governs the bending rigidity. Red blood cells have 
a bending rigidity Eb ~ lO/c^T, an area stretching mod- 
ulus K ^ EHq ^ 10/iN/m [28], and sizes Rq ^ 4 /im, 
which results in a smaller dimensionless bending modu- 
lus Eb ^ 0.0005 and a much smaller classical buckling 
threshold \pcb\ ~ 0.2 Pa in accordance with the fact that 
red blood cell shapes are buckled at ambient conditions. 

From our above result Pc ~ 0.23 Pcb for thin shells with 
Eb ^ 0.001, we expect critical values \pc\, which are 
smaller by a factor of at least 5 compared to the classical 
buckling pressure \pcb\ for all of these capsules. 



VI. ANALYSIS OF SIMPLY BUCKLED SHAPES 
AND BENDING MODULUS 

Smaller bending resistances allow sharper bends in 
buckled configurations. Hence, the minimal radius of 
curvature l/^^edge = niins(l//^s) (in 5-direction along the 
contour), which occurs close to the edge of the indenta- 
tion of a simply buckled shape 2 (see figure 12), should 



depend sensitively on the reduced bending resistance Eb 
and represent an adequate observable to infer the reduced 
bending modulus Eb- 

Figure [TT] shows a double logarithmic plot of the mini- 
mal radius of curvature as a function of the reduced bend- 



14 



(l/^smax)/^o 
0.5 




Figure 11. (Color online) Minimal radius of curvature l/^^edge 
plotted against the reduced bending resistance Eb- Dots: 
Computed shapes for V = 0.8 Vb. Line: Power law fit for 
Eb < 0.004. 




Figure 12. (Color online) Geometry of an axisymmetric dim- 
ple in a simply buckled shape 2. 



ing modulus. It was obtained from a series of simply 
buckled shapes with fixed volume V = O.SVq. For bend- 
ing moduli Eb < 0.004, a power law can be fitted to the 
data: 



1 



D ^0.2356 



(37) 



The power law ( 37 ) can be confirmed by a scaling argu- 
ment, where we balance bending and stretching energies. 
We consider a small dimple with radius r <C i^o and 
depth h as depicted in figure [12] The elastic energy is 
mainly located at the edge of this dimple, which spans 
the width d <C r. 

Since the dimple is assumed to be mirror-inverted to 
its original shape, its radius and depth are determined 
by the angle a (see figure 12). To leading order, they 
are given hy r ^ Roa and h ^ Rqo^ respectively. By 
that, the volume change with respect to the unbuckled 
configuration is approximately AV h r'^ Rq . 

The energies of bending and stretching are of the order 



and 



EHo^rd 



(38) 



respectively [13j, where ( is the typical radial displace- 
ment of the membrane near the edge of the dimple. 
Because the direction of the meridian changes about a 
within the width we have ( ~ ad ~ d (AV/Rq)^^^. 
The total elastic energy for a given volume reduction A V 



therefore takes the form 
F - EHo- 



R. 



7/2 



Eb- 



'AFi/2 



dR, 



3/2 



(39) 



Minimizing the total elastic energy with respect to d we 
find the equilibrium width of the dimple edge. 



EHc 



Rt 



1/4 



^0 E^B^^ 



(40) 



which can also be written in the form d ~ \/Hq Rq as 
in ref. [13j. Confining the directional change a of the 
meridian to a width d of the edge of the indentation (see 
figure 12) results in an edge curvature 



^edge 



a 
d 



AVi/4 



-1/4 



R. 



7/4 



(41) 



When the critical buckling pressure \pc\ is small (com- 
pared to the pressure unit defined by EHq/Rq), the un- 
buckled region outside the dimple remains roughly spher- 
ical with a radius close to the original radius Rq, as can 



be seen from eq. 31 Therefore 



AV = Vb- V = (l-v)Vb 



(42) 



holds to a good approximation, where v = V/Vq is a 



reduced volume. Using this in (41 ), we find a scaling law 



^edge ~ 5 E^ , 



(43) 



which is of the same form as the above fit (37) with an 



exponent 1/4 matching the fit result 0.2356 from (37) 
quite well. 

It is evident from the assumption d <C r, which cor- 
responds to a sharp edge of the indentation, that the 
scaling law holds only for sufficiently small bending resis- 
tances, Eb < 0.004 in this case. The assumption r <C i^o 
implies that the scaling law holds for sufficiently small 
volume changes. Indeed, analyzing the scaling behavior 
for several capsule volumes using (43), we find that the 



fitted exponent matches the theoretical value 1/4 very 
well for V = 0.9 Vb or 0.8 Vo, but starts to deviate for 
V = 0.67 Vo or 0.51^, see table ITT) We can also deter- 



mine that the numerical prefactor in (43) is of the order 



of unity and only weakly volume dependent. In contrast 
to these findings for simply buckled shapes 2 without 
contact of opposite sides, we observe that /^edge is nearly 
independent of Eb for buckled conformations of branch 
2c with opposite sides in contact. 

The results of the fits presented in table [IT] could be 
used to quantitatively analyze experimental shapes of 
simply buckled elastic capsules without opposite sides in 
contact as shown, for example, in Refs. [8l[T0l[T2l[T8], pro- 
vided the radius of curvature at the edge of the dimple 
can be measured accurately. We note that the calculated 
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reduced Volume v 


prefactor c 


exponent b 


0.9 


0.73 


0.246 


0.8 


0.63 


0.236 


0.75 


0.59 


0.229 


0.67 


0.54 


0.221 


0.5 


0.50 


0.217 



Table II. Fit parameters for different capsule volumes. The 
fit model for the radius of curvature at^the dimple edge is 

l//^edge = C(l 



-1/4 



Ro , see eq. ( 43 ) . For small volume 



changes, the theoretical exponent 1/4 agrees best. 



determine the critical pressure Pc- 

For buckled shapes, the maximal curvature /^edge at 
the edge of inward buckled dimples was found to de- 
pend depend strongly on the ratio of bending resistance 
to surface Young modulus, with smaller ratios leading 
to sharper bends. A power law /^edge ^ {EHq/ EbY^^ 
was found for sufficiently small bending resistances and 
sufficiently small volume changes. This relation may be 
used to analyze experimental shapes of buckled elastic 
capsules and extract the bending modulus of the capsule 
membrane from the capsule shape. 



shapes as shown in figure [TT] are in qualitative agree- 
ment with some of the experimentally observed shapes 
[HI Uni [121 [IHl- From a measurement of the curvature 
at the edge of the buckling dimple A^edge the dimension- 
less bending modulus Eb = Eb/{Rq E Hq) can be de- 
termined using (41 ) and the numerical prefactor from ta- 



ble |II[ In combination with an independent measurement 
of Young's modulus E^ for example, via a measurement 
of the classical buckling pressure Pcb, this type of shape 
analysis provides a method to obtain the bending mod- 
ulus of a capsule, which is hard to measure otherwise. 



VII. CONCLUSIONS 

We applied nonlinear shell theory to the problem of 
axisymmetric deformations of an initially spherical cap- 
sule. The elastic properties of the capsule membrane 
were modeled with a quadratic strain-energy function. 
This approach to Hooke's law allowed us to use the meth- 
ods of force and moment equilibrium and a least-energy 
principle simultaneously. 

Bifurcation diagrams for reduced capsule volume and 
negative pressure were presented. The least-energy prin- 
ciple gave information about the preferred configurations 
and allowed to obtain a complete picture of the transi- 
tion into the fully buckled state. If the capsule volume 
is controlled, simply buckled configurations turned out 
to be energetically preferable below a certain critical vol- 
ume Vc, but an energy barrier must be overcome in order 
to leave the trivial solution branch. The transition at 
Vc is thus discontinuous. If the internal negative pres- 
sure is controlled, spherical shapes become mechanically 
unstable at the classical buckling pressure Pcb- However, 
buckled shapes with opposite sides in contact become en- 
ergetically favorable at a much lower negative pressure Pc^ 
i.e., \pc\ < \Pcb\' Also at controlled negative pressure, the 
transition at Pc is discontinuous and involves an energy 
barrier. With the methods presented here, also configu- 
rations with opposite sides in contact could be computed 
and incorporated in the bifurcation diagrams; fully buck- 
led configurations with opposite sides in contact (branch 
2c in the bifurcation diagrams) actually have the lowest 
energies at small volumes or low negative pressures and 
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Appendix A: Calculus of Variations 

In this appendix we derive the first variation and the 
resulting Euler-Lagrange equations for the enthalpy func- 
tional 



G 



(27TroWs—p7Tr'^XsSimjj^ dsQ. (Al) 



The integrand has to be regarded as a function of sq. In 
functions like r{s), a change of variables from s to sq can 
be performed by the function s{so) introduced in ([5|. 

Two of the functions r(5o), ^(^o) and ^^(so) determine 
the capsule configuration completely. As the integrand 
of (14) and the geometrical relations contain mainly r 
and ip, we choose r{so) and tp{so) as the two basic fields. 
Variations Sr{so) and SiIj{so) have to fulfill the boundary 
conditions 

Sr{0) = Sr{Lo) = and SipiO) = SipiLo) = 0. (A2) 
The first variation SG of G[r, t/j] is obtained as 

SG = / {27r To Sws — pTT S{r'^ Xg sinilj)) dsQ. (A3) 
Jo 

The variation 6ws of the surface energy density intro- 
duces tensions and bending moments with the help of 
the constitutive equations (10), 



dws ^ , dws ^ , dws ^ ^ , dws . ^ 



dcA 



OK, 



= Ts Scs + Xs T(i, Se^f, + rus SKg + A^ 5K^. 

(A4) 

Now, the variations of the strains must be expressed in 
terms of (5r, 5ip and its derivatives with the help of strain 
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definitions and geometrical relations, 



^ cos?/^(so) 
es = 1 

sin 

= - Ho 



SKs 



5r' 

— Sr 
ro 

ro 



As tan ip 5iIj 



For a stationary shape, 5G = for arbitrary variations 5ip 
and 5r^ and the terms in curly braces have to vanish. This 
gives the Euler-Lagrange equations describing stationary 
states. Rearranging the term next to 5iIj by a change of 
variables ds = Xg d^o, we obtain 



(A5) 

In a similar fashion, the variation of the second term can 
be calculated as 



with q 



tan?/; 



(5(r^ \s sin ?/^) = 5{r' tan ?/^) 

= 2r Xg sin 5r ^ tan ijj Sr' 



1 d(r rris) 

r ds 

1 r 

' T' 

2 cosip 



(A8) 
(A9) 



(A6) 



Inserting everything into ( [A3| ), and sorting according to which are eqs. ([T5| and (frfl) in the main text. With this 



Sr, Sr', Sip, Sil)' yields a rather long expression. Using 
integration by parts to transform Sr' into Sr and Sip' 
into Sip (note that the boundary terms vanish) results in 
the following result for the first variation of G, 



definition of the transverse shear tension q, the term next 
to Sr can be simplified to 



SG = j d^o ^(^r|27r As — 27rpr As sin?/^ 

d f rTs \ d , 9 ,x1 

—-4 T]^^P-^y^ t^nip)\ 

dso^co^ipJ dso J 



27r- 



1 d(rrs 
r ds 



cos tp 



(AlO) 



Sip^^27rrTs As tan ip -\- 27r Xs rricj) cos ip 



cosip 



27T- 



which gives eq. (16) in the main text. 
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